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ABSTRACT 


This report presents a detailed description of a computer 
program to calculate cumulative damage of notched structural members 
subjected to arbitrary spectra. The Sequence Accountable Fatigue 
Analysis computer program develops its sequence sensitivity by 
tracking residual stresses local to a notch throughout the spectrum 
of loads. Residual stress relaxation analysis is included to increase 
the generality of the results. An example spectrum and resulting 


cumulative damage analysis are illustrated. 
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SYMBOLS 


Residual stress 

Maximum local stress level 

Minimum local stress level 

Yield stress 

Equilibrium component of the residual stress 

Local strain, total 

Elastic component of total strain 

Plastic component of total strain 

Strain intercept at one reversal on a log 5 198 life curve 

Slope of the log Ee, 108 life curve 

Maximum applied stress level 

Minimum applied stress level 

Mean applied stress level, (Snaxt Smind/2 

Elastic stress concentration Factor 

Damage 

Integer describing the level number 

Equilibrium period, number of cycles for the local 
Stresses to return approximately to the equilibrium 
conditions following an overload 

Residual stress relaxation constant 

Relaxation function exponents 

Number of cycles 


Modulus of elasticity 


Number of cycles of life at a given stress or 
Strain cycling level 
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SECTION I 


INTRODUCTION 


Cumulative damage analyses based upon the local stress-strain 
behavior at a notch appear to be reasonably successful in anticipating 
trends in fatigue life behavior of notched specimens subjected to 
spectrum loading (1-6), The type of ee that usually occurs is 
that peak tensile loads tend to increase the fatigue life and peak 
compressive loads tend to decrease the life of notched structures 
compared to structures experiencing load spectra not having those 
peaks (5,6) | Local behavior analyses, such as those developed by 
Smith (7) and Neuber (8) help to explain this phenomenon as being 
a result of the tensile peak load creating a compressive residual 
stress at the notch and, conversely, the compressive peak creating 


a tensile residual stress. The change in life occurs because the 


-residual stress state modifies the subsequent damage accumulation 


rates. 

The Sequence Accountable Fatigue Analysis computer program was 
developed to incorporate the local stress-strain approach with a 
recent residual stress relaxation analysis (6) order to improve 
the sequence sensitivity of cumulative damage analysis. This 
technical memorandum presents the details of the resultant computer 
program and an example of its use. The correlation of predictions made 
with this analysis to actual results of tests epper leneine spectrum 


(9) 
? 


loading is presented by Potter and Potter, Gallagher, and 


Stalnaker (10). 


SECTION II 


PROGRAM OUTLINE 


The Sequence Accountable Fatigue Analysis traces the stress-strain 
behavior local to a notch throughout an applied load spectrum and 
calculates the damage based on the local experience. The computer 
program is divided generally into the four parts or modules outlined 
in Fig. l. 

The basic input data for the material, specimen geometry, fatigue 
behavior qualities and spectrum, are developed in Module I. The infor- 
mation required in Module I is discussed further in Section III. Module 
II takes the input information and determines the local stress~strain 
behavior. Module III references the Range Pair Counting Method Sub- 
routine to cycle count the local stress spectrum. Module IV determines 
the damage in the local stress-strain spectrum. 

The basic analyses used in Modules II, III and IV are presented 
below. 

Module II - Local Stress-Strain Behavior - The analysis used during 
the determination of the local stress behavior during the spectrum 


(7) 


of loading is a combination of analyses developed by Smith 


(8) (6) 


Neuber and Potter Smith's simple analysis indicated that 
the residual stress could be approximated by assuming that the initial 
stress-strain behavior was elastic upon unloading following plastic 


flow. Thus, the residual stress could be calculated knowing the 


maximum local stress and the maximum applied stress as in Eq. 1 and 


in Fig. 2. 


Ova. Ones. = ReSaase. (1) 
res, Max t max; 


The Sequence Accountable Fatigue Analysis computer program currently 
incorporates elastic-perfectly plastic stress-strain behavior. There- 


fore, Dao is equal to the yield stress. For the cycles immediately 
i 


following the peak stress, the residual stress determined in Eq. l 
modifies the elastic solution as shown in Eqs. 2 and 3 (provided that 


the following maximum applied stress is less than Dae and that there 


is no change in the residual stress due to a minimum applied stress 


causing reversed yielding). 


oF O + K.-S 2 
max. res._4 t max; (2) 


Onin, ress_4 - Ke Smin, (3) 


The analysis developed by Neuber (8) has been extended to cyclic 


(3) 


loading by Wetzel (2) and Wetzel, Morrow and Topper and used by 


(1,4-6 


) 
many others primarily to determine local stress-strain 
behavior. It is used in this program only to calculate plastic strains 
occurring when the residual stress undergoes a step change. The 


plastic strain calculation routine is accessed only when the Oe: 
i 


ord a terms in Eqs. 2 and 3 exceed tensile or compressive yield 
A. 
stress levels, respectively. Figure 3 illustrates the calculation 


of the plastic strain. 


(2) 


The local stress~strain behavior, according to Wetzel is 
related to the applied load by Eq. 4 


)*/E (4) 


ote = (KS 
t max 


The plastic strain can be found by subtracting the elastic com- 


ponent from the total strain. 


Zz 
[pe 7b =e (Kt Smax) /E*Omax 7 Smax/E 
Therefore, the plastic strain associated with Da is given in Eq. 5. 
i 
ae 2 - 
Dy = (KtSmax,) [EOy 5 Oys/E (5) 


If a residual stress existed prior to this plastic strain 
excursion, the plastic strain associated with that prior excursion 


is subtracted from Eq. 5 as shown in Eq. 6. 


yi 2 
“pay = (Ki Sax, ) [Edys ~ (yg [res;_y) [Edy 5 (6) 
A similar calculation is made for plastic strains occurring during 


the minimum stress peak. 


(6) 


In the analysis developed by Potter the residual stress 


cyclically rehaxes toward zero or an equilibrium residual stress as 
shown in Fig. 4 according to Eq. /. 


0 )exp(N/Ngp, 1n(0.1)) (7) 


eSN=1,2,... — ‘Sresya1 7 *resxQ 
The Nop term, the Equilibrium Period, is dependent upon the 


applied stress and the Residual Stress Relaxation Constant. 


El E2 


The Residual Stress Relaxation Constant, Cl, has not yet been 
experimentally defined but should be a constant for a material. 
Module III - Cycle Counting Method 
After the local stress and plastic strain behavior is calculated, 
the local stress spectrum is Range Pair Counted using a computer program 


developed by Tigchier. +? 


Module IV — Damage Calculation 

Damage is calculated separately for the plastic strain excursions 
and the elastic stress spectrum. The damage is determined from the 
conventional D =n coteuiecien: Damage from each of the plastic strain 


12 
cycles is determined from the Coffin-Manson (12) form 


Dy = Le /Ng, = 1-/(ep,/egr) 


Damage from the elastic stress cycles is determined in a similar manner. 
The maximum and minimum local stress levels are sequentially compared to 
unnotched S-N data in a Modified Goodman Diagram format. Damage is 

summed, and failure of the coupon is defined as the event occurring when 


the summed damage equals unity. 


SECTION III 


INPUT DATA REQUIREMENTS 


In general, each spectrum analyzed will require slightly different 
programming in order to get the load history into a usable format for 
the core program. The basic program requires a certain family of in- 
formation before any analytical predictions can be made. Appendix I 


contains a program listing for the Sequence Accountable Fatigue Analysis. 


The subroutine CORE which accesses the subroutines having to do with RPCM, 


the Range Pair Counting Method, contains the basic analysis. Subroutine 
SAL reads the data input and then references subroutine CORE. The 
subroutine SAL shown is one in which a block of eveles te repeated with 
optional cycles. A list of the input data cards and the resulting 
analysis is given in Appendix II. 

The specific data requirements are given below. 

1. Stress-Strain Behavior - The stress-strain behavior is presumed 
to be elastic-perfectly plastic with the tensile yield stress being 
equal to the compressive yield stress. The yield stress value used is 
an average of the monotonic behavior generally being above the 0.22% 
yield value and below the engineering ultimate strength. 

2. Residual Stress Relaxation - The residual stress relaxation 
behavior of Eq. 7 and 8 is characterized by Cl, the Residual Stress 
Relaxation Constant and El and E2, the relaxation equation exponents. 
The Residual Stress Relaxation Constant, Cl, has not yet been adequately 


determined. It should be a material property if the relaxation function 


is correct and must be assumed. A reasonably accurate estimate of the 
Residual Stress Relaxation Constant for aluminum material falls in the 
range of 5-20 x 10° (cycles) (Ksi)*. Further experimentation on the 
part of the analyst should develop a Cl usable for his set of conditions 
until actual measurement of residual stress relaxation behavior defines 
the relaxation function and constants. At present El and E2 are considered 
to be equal to 1.0. 

3. Specimen Geometry - The elastic K, value (if available) is 
entered into the analysis. If that value is not available then an 
estimate from some other method may be used. In certain cases, a value may 
be determined from a constant amplitude fatigue test of a similar struc- 
ture by fitting several values of Kt to the analysis and determining 
the best correlation as is done with the Ke solution. Once a stress 
concentration factor, K,, is determined for a specimen, that value is 
not changed from test-to-test of the same coupon donfiguration. 

4. Load Multiplier - Different spectra are presented for analysis 
in different manners. Some data are presanted in percent of maximum 
stress, others in terms of nominal stress, and others in terms of bending 
moment. The value of the load multiplier defines the nominal stress 


history. 


5. Cumulative Damage Analysis - The damage from the range-paired 
n 


elastic stress spectrum is determined by calculating a simple N value 


for each level and accumulating the total. The Mes value is determined 


from unnotched coupon S-N data in the Modified Goodman Diagram format. 


The program requires the input of four second order equations describing 
the maximum and minimum stress levels at lives of 10", 16", 10° and 10/ 
cycles. The coefficients of the equations are derived by least square 
fitting the S-N data presented in the form of Eq. 9. 


Ss =A(I)S_, 2 4B(I)S_. + C(I) (9) 


max in 
The A, B, and C coefficients for several typical materials are presented 
in Appendix IV. The S-N data shown was derived from various sources 


(13) 


but usually from the MIL-HDBK-5A The C coefficients correspond to 


the maximum stress level at zero to maximum applied stress conditions 


on the unnotched coupons. 
The damage from the plastic strain cycles is determined using the 


Coffin-Manson relation to calculate the Ny value. The conventional 


e 


i 
plastic strain intercept at one reversal and the a life slope values 


are used in the analysis. Specific measured values from the literature 
are used when available and typical values when they are not available. 
6. Analysis or Test Spectrum - The last information needed is the order 


and magnitude of application of the spectrum used in the test. 


SECTION IV 


OUTPUT OPTIONS 


The Computer Program prints the following output in the process 
of the analysis. 

1. Maximum and minimum applied stress and local stress response 
through the spectrum. Also printed out is the residual stress, equilib- 
rium stress, applied cycles, and the equilibrium period. 

2. The elastic local stress history as input into the Range Pair 
Subroutine and the resulting Range Paired spectrum. 

3. The plastic strain occurrence during the spectrum and the damage 
associated with each strain reversal. 

4. The accumulated ee associated with the plastic strains. 

5. The Range Paired elastic stress spectrum and the damage 
associated with each level. 

6. The accumulated damage associated with the current block of 
loading including the plastic strain damage and the total damage shee the 
initiation of cycling. 

At the option of the analyst, he can print out all the above items 
or only two. The IPRINT value controls what data is printed. 


If IPRINT = 1, all six items are printed for each flight or block. 


If IPRINT = 2, all items except 2. above are printed. 


If IPRINT = 3, only items 4. and 6. above are printed. 
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APPENDIX III 


LIST OF COMPUTER PROGRAM SYMBOLS AND DEFINITIONS 


A Coefficient of the x? term in the equation of a line on 


a constant life fatigue diagram where minimum stress is x 


and maximum stress is y. (R= Ax” + Bx + C - y) 
AA An assigned value of +1. or -1l. 
AAA A stress used in the calculation of plastic strain. 


ABDIF The absolute value of DIF. 

ABM The absolute value of ASMAX or of ASMIN, as assigned. 

ABMAX The absolute value of ASMAX. 

ABMEAN The absolute value of ASMEAN. 

ABMIN The absolute value of ASMIN. 

ABR4 The absolute value of R(4). 

ABR7 The absolute value of R(7). 

ABS The name of a routine calling for the absolute value 
of a quantity. 

AKT Stress concentration factor, K.. 

ASMAX The product (ART) (STMAX) 

ASMEAN The quantity (ASMAX + ASMIN) /2 

ASMIN The product (AKT) (STMIN) 

AVSGMN Average value of SIGMIN over an interval. 


AVSGMX Average value of SIGMAX over an interval. 
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B Coefficient of the x term. (See A.) 
BBB A stress used in the calculation of plastic strain. 
C The constant. (See A.) 


COFMAN Inverse of the Coffin-Manson slope. 
CYCINT The number of cycles in an interval. 
CYCLES The calculated number of cycles expected to be indicated 
on a constant life fatigue diagram for the applied combination 


of maximum and minimum stress. 


Cl The Residual Stress Relaxation Constant (See ENEP.) 
DAM Damage. 

DECK Decimal or real value of integer K after conversion. 
DEL2 A portion of a least-squares-method solution. 

DIF The difference between residual stress and equilibrium 


residual stress. (RES(1I) -— EQRES) 

DO2 A portion of a least-squares-methnod solution. 

DUMMY A variable used in the calculation of the number of cycles 
to be considered as an interval for relaxation determination. 

ELMOD The elastic modulus. 

EN The number of cycles from the beginning of the relaxation 
process to the end of the current interval. 

ENEP The number of cycles required for overload residual stress 
effect to return to within one-tenth of its original difference 
from equilibrium conditions. 


( Ng, = CL/ (ABM) El (apMEAN)=* ) 
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ENN 


ENNCYC 
EPSD 
EQRES 
EX 


EXP 


EXPO 


IBLOCK 


IFIX 


IN 


IPRINT 
TRAIN 


IRPCM 


The number of applied cycles at a load level. 

The ratio of the number of applied cycles to the number 
of cycles to failure. ( ENN/CYCLES ) 

LCF strain intercept. 

Equilibrium residual stress. 

An exponential function depicting the relaxation of 
residual stress. 

The name of a routine calling for the exponential value 
of a quantity. 

An exponent. The power of 10 which indicates the number 


of cycles to failure. 
Residual stress Relaxation Exponents. 


The name of a routine calling for integer-to-real conversion. 
A variable subscript. 

The iaeued fein aimed of a block the blocks being 

numbered consecutively from 1 to NBLOCK. 

The name of a routine calling for real~to~-integer conversion. 
The number of steps input to the range pair counting 
subroutine. 

Value controlling the WRITE statements. 

A counter. 

Value controlling entry into the range pair counting 


subroutine. 
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ISTEP 


ITYPE 


KPMAX 


LMN 


NBLOCK 


NDECK 


NFLAG 


NFLAG2 


The identifying step number, the steps being numbered 
from 1 to NLEVEL. 

The identifying type number, the types being numbered 
from 1 to NTYPE. 

A variable subscript. 

Value of +1 or 0, as assigned for branch determination. 
Value of -l or 0, as assigned for branch determination. 
An index variable. 

An index variable. 

An index variable. 

An index variable. 

An index variable. 

The number of steps output from the range pair counting 
subroutine. 

An index variable. 

An index variable. 

An index variable. 

An index variable with values of N=4-7 indicating the 
power of 10, and thus identifying a particular life cycle 
curve. 

The total number of times to execute a block of loads. 
The number of data decks to be run sequentially. 

An integer used as a counter. 


An integer used as a counter. 
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NLEVEL 


The total 


number of steps, or levels, of loads in a block. 


A subscripted variable used to indicate which types of 


experienced in which blocks. 


number of different types. 


Residue term in damage calculation. 


The number of cycles for a level after exitting the range 


of damage for a flight. 


Accumulated total of applied cycles. (Summation of ENN). 


Accumulated cycle ratio. (Summation of ENN/CYCLES). 


of R(N), N=4,7. 
of NR(N), N=4,7. 
of R(N)*, N=4,7. 
of NR(N) 7, N=4,7. 
of R(N)”, Ne4,7. 


of R(N)*, N=4,7. 


Identification of the source of the SN data. 


NN 

loads are 
NTYPE The total 
PLSTRA Plastic strain. 
R 
RES Residual stress. 
RNCYC 

pair counting subroutine. 
S IGMAX Maximum stress. 
SIGMIN Minimum stress. 
STMAX Maximum applied stress. 
STMIN Minimum applied stress. 
SUMDEL Summation 
SUMENN 
SUMNC 
SUMR Summation 
5S UMRN Summation 
SUMR2 Summation 
SUMR2N Summation 
SUMR3 Summation 
SUMR4 Summation 
TITLE] , TITLE2 
TLL Tensile load limit. 
TM1 , TM2 Material type. 
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TTYS One«fifth of tensile yield stress. 
TYS Tensile yield stress. 


T1,T2,T3,T4,T5,T6,1T7,T8 Test identifying information. 


x Variable equivalent to SIGMIN. 


y Variable equivalent to SIGMAX. 
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MATERIAL 
2024-T4 


2219-T851 


7075~T6 


RQC-100 


Man~Ten 


4340 Steel 


Ti-6-4 


CO CO O 


APPENDIX IV 


FATIGUE LIFE INPUT DATA FOR SEVERAL MATERIALS 


YIELD 


STRESS 
58. 


55. 


72. 


125. 


55. 


160. 


158. 


STRAIN 
INTERCEPT 


0.4 @ 


0.4@ 


0.4 © 


0. 42) 


1.1 


0.4@ 


0.4@ 


INVERSE 
OF SLOPE 


-1.8362) 


-1.8362) 


-1.836D) 


1.4958) 


-1.66%)) 


1,83) 


-1. 8362) 


Division, Courtesy of Mr. H.R. Jaekel. 
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LIFE, 


104 


SIS UT SD “SO Ul & SON UT “ON UT 


SAM 


S-N LIFE COEFFICIENTS 


A(T) 


~.0020 
-.0032 
~.0035 
~.0042 


-.0022 
~.0018 


-,0015 
~.0024 


-.0020 
-.0022 
-.0014 
~.0013 


-.0002 
~.0007 
-.0005 
-.0005 


-.0009 
-. 0006 
-.0000 

0001 


B(T) 


2091 
- 4366 
6207 
- 7003 


2204 
- 3320 
4628 
6420 


2801 
5154 
6141 
6838 


- 2136 
2927 
- 3669 
- 4376 


2257 
» 3520 
- 4669 
5678 


2567 
5248 
~5557 
25957 


. 2368 
. 4640 
4650 
4752 


C(I) 


62.6@) 


51.4 
42.2 
36.1 


55.84) 


48.3 
39.7 
31.7 


71.78) 


56.3 
44.6 
38.1 


98.) 


88.5 
79.1 
70.3 


63.50) 


53.1 
43.7 
35.4 


162.4) 


126.9 
113.5 
108.5 


154.23) 


110.3 


88.9 
84.2 


Data from Endo, T., and Morrow, J., NAEC-ASL-1105, Naval Air Engineering 
Center, Philadelphia, PA, June 1966. 

Data not available - Source 1 data considered typical. 

Derived from Metallic Materials and Elements for Aerospace Vehicles Struc- 
tures, MIL-HDBK-5A, Dept. of Defense, Washington, D.C., February 1966. 

Information supplied from Rockwell International. 

Information supplied by Society of Automotive Engineers Cumulative Damage 
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